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ABSTRACT 

We compute the real-space cluster-galaxy cross-correlation £ cg (r) using the 
ROSAT-ESO Flux Limited X-ray (REFLEX) cluster survey, a group catalogue con¬ 
structed from the final version of the 2dFGRS, and galaxies extracted from 2MASS 
and APM surveys. This first detailed calculation of the cross-correlation for X-ray clus¬ 
ters and groups, is consistent with previous works and shows that £ cg (r) can not be 
described by a single power law. We analyse the clustering dependence on the cluster 
X-ray luminosity Lx and virial mass M V i r thresholds as well as on the galaxy limiting 
magnitude. We also make a comparison of our results with those obtained for the 
halo-mass cross-correlation function in a ACDM N-body simulation to infer the scale 
dependence of galaxy bias around clusters. Our results indicate that the distribution 
of galaxies shows a significant anti-bias at highly non-linear small cluster-centric dis¬ 
tances ( b cg {r ) ~ 0.7), irrespective of the group/cluster virial mass or X-ray luminosity 
and galaxy characteristics, which show that a generic process controls the efficiency 
of galaxy formation and evolution in high density regions. On larger scales b cg (r) rises 
to a nearly constant value of order unity, the transition occuring at approximately 
2 h -1 Mpc for 2dF groups and 5 h -1 Mpc for REFLEX clusters. 
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1 INTRODUCTION 

Statistical analysis of clusters and groups of galaxies are 
among the most powerful tools for the study of large scale 
structure. In the hierarchical clustering scenario, these ob¬ 
jects form by the gravitational amplification of primordial 
density fluctuations, so that they can be used to characterize 
the spatial distribution of the peaks of the density field. Be¬ 
sides, as they are the largest gravitationally bound objects 
and span a wide range of masses, they can be used as labo¬ 
ratories to determine the role of the different astrophysical 
processes that govern galaxy formation. 

The first statistical analysis on this type of objects 
where based on cluster catalogues constructed by visual 
identification (Abell, 1958; Abell, Olowin & Corwin, 1989). 
The cluster-cluster correlation function £ cc (r) has been the 
preferred tool for such analysis, providing a simple and con¬ 
venient statistical tool to characterize their spatial distribu- 
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tion and place constraints on cosmological models (Bahcall 
& Soneira 1983; Klypin & Kopylov 1983; Peacock & West 
1992; Postman, Huchra & Geller 1992). Nonetheless, clear 
evidence of inhomogeneities and line-of-sight projection ef¬ 
fects was found (Lucey 1983; Sutherland 1988), which where 
responsible for part of the observed amplitude of the angu¬ 
lar correlation. The situation improved with digitized clus¬ 
ter surveys like the Edinburgh/Durham Cluster Catalogue 
(Nichol et al. 1992) and the APM survey (Dalton et al 1992, 
1994; Croft et al 1997). Although these where more homoge¬ 
neous, systematics due to projection effects are inherent to 
cluster identification based on angular positions in the sky 
(Padilla & Larnbas, 2003). 

In recent years, X-ray information has been used to con¬ 
struct new cluster samples. This method of identification has 
several advantages with respect to the use of optical data. 
The X-ray emission is a strong signature of a gravitational 
potential well since there is a strong relation between the X- 
ray luminosity Lx and mass (Reipricli & Bohringer 2002). 
Besides the emissivity of thermal Bremsstrahlung radiation 
is proportional to the square of the electron density thus 


© 0000 RAS 



360 



0 


Figure 1 . a) Aitoff projection of the distribution, in galactic coordinates, of 2MASS galaxies with K s < 12.5 (dots) together with the 
REFLEX cluster sample (red circles). The angular mask described in section 2.2 has been applied to both samples, b) Aitoff projection 
of the distribution, in equatorial coordinates, of 10% of APM galaxies with bj < 18.5 (black dots) and all 2dFGGC (blue dots). 


it provides clear evidence of a mass concentration which is 
not strongly affected by projection effects. Furthermore, the 
X-ray emission from the clusters is concentrated towards 
the dense central core giving an improved angular resolution 
compared to the one obtained using galaxy concentrations, 
also helping to reduce the possibility of projection effects. 
The ROSAT-ESO Flux Limited X-ray (REFLEX) cluster 
survey (Bohringer et al. 2004), comprising 447 objects, is 
the largest statistically complete X-ray selected sample. It 
has been previously used to analyse the large scale struc¬ 
ture through its power spectrum (Schuecker et al. 2001) 
and its correlation function (Collins et ah, 2000). Although 
Schuecker et al. (2001) found a systematic increase of the 
amplitude of the power spectrum with limiting X-ray lu¬ 
minosity, Collins et al. (2001) results lack such a trend in 
the correlation function, characterized by an approximate 
power-law with correlation length r o = 18.8 ± 0.9 h -1 Mpc. 

On the other hand, since the advent of large galaxy 
redshift surveys a reliable identification of groups of galax¬ 
ies has been possible using not only angular positions, but 
also the complete redshift information, avoiding in this way 
a strong contamination by projection effects. The first algo¬ 
rithms designed for this task where developed by Huchra & 
Geller (1982) and Nolthenius & White (1987). These meth¬ 
ods where based on the friends-of-friends algorithm used 


in numerical simulations with a redshift dependent link¬ 
ing length parameter. A variety of group catalogues where 
constructed using this algorithms: Merchan et al. (2000) 
used the Updated Zwicky Catalogue (Falco et al. 1999), 
Giuricin et al. (2000) constructed a group catalogue from 
the Nearby Optical Galaxy Sample, Tucker et al. (2000) ex¬ 
tracted a group sample from the Las Campanas Redshift 
Survey (Shectman et al. 1996), and Ramella et al. (2002) 
used the information of the Updated Zwicky catalogue and 
the Southern Sky Redshift Survey (da Costa et al. 1998). 
Recently, the 100K public data release of the 2dF Galaxy 
Redshift Survey (2dFGRS, Colies et al. 2001) has been used 
by Merchan & Zandivarez (2002) to construct one of the 
largest group catalogues using a slight modification of the 
Huchra & Geller (1982) algorithm. The large scale redshift 
space distribution of this sample has been studied by Zan¬ 
divarez, Merchan & Padilla (2003). They showed that the 
redshift space correlation function has a power law behavior 
with a correlation length of so = 8.9 ± 0.3 /i _1 Mpc. 

The cluster and group catalogues identified by X-ray in¬ 
formation or three dimensional percolation algorithms, rep¬ 
resent different types of objects. In order to be able to con¬ 
front the results obtained for these samples with the pre¬ 
dictions of theoretical models, a complete understanding of 
the differences in the statistical properties of the samples is 
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crucial. One of the most important possible comparisons is 
the analysis of the clustering of galaxies, selected in differ¬ 
ent wave-bands, around these objects and how it depends 
on the properties like X-ray luminosity or virial mass. The 
best statistical tool for such analysis is the two-point cluster- 
galaxy cross-correlation function £ C9 (r), which is a measure 
of the mean radial density profile around clusters. The first 
measurement of this correlation function was performed by 
Seldner & Peebles (1977), who studied the angular cross¬ 
correlation of Abell clusters and Lick galaxies. They found 
that £ cg (r) was very large and positive out to scales as large 
as r ~ 100 /i _1 Mpc but they where contradicted by Lilje 
& Efstathiou (1988) who did not find any strong evidence 
of clustering beyond r >= 20 /i _1 Mpc arguing that part of 
the signal detected by Seldner & Peebles (1977) was a conse¬ 
quence of artificial surface density gradients in the Lick cata¬ 
logue. At smaller scales, Lilje & Efstathiou (1988) found that 
£,cg(r) was well fitted by a power law with ro = 8.8 fe _1 Mpc. 
Later works (Dalton 1992; Mo, Peacock & Xia 1993; Moore 
et al. 1994) found similar results based on direct measure¬ 
ments from redshift surveys of galaxies and clusters. They 
found that £ cg (r) has a similar shape to £ cc (r) and £ gg (r) 
but with an amplitude approximately equal to their geo¬ 
metric mean. Croft et al. (1999) estimated £ cg (r) in real and 
redshift-space from the APM galaxy and cluster surveys and 
found that its shape can not be described by a single power 
law and that its amplitude is almost independent of cluster 
richness and the limiting magnitude of galaxies. 

In this paper we compute £ cg (r) using a sample of 
groups derived from the final version of the 2dFGRS and the 
REFLEX cluster survey, and the galaxies from the 2MASS 
(Skrutskie et al. 1997) and APM (Maddox et al. 1990) sur¬ 
veys. The outline of the paper is as follows, in section §2 we 
describe the different data sets analysed in this work. In §3 
we describe the method used to obtain the 3D real space 
cluster-galaxy cross-correlation function from the projected 
cross-correlation function. In §4 we test the dependence of 
our results on the magnitude limit in the galaxy catalogue, 
the X-ray luminosity Lx in the REFLEX catalogue and the 
virial mass M vlr in the 2dF Galaxy Group Catalogue. In §5 
we show a comparison with the results of an N-body simu¬ 
lation of the ACDM cosmological model. Finally, in §6 we 
present a short discussion and our main conclusions. 


2 THE DATA 

2.1 Galaxy catalogues 

We used the information from the Two Micron All Sky Sur¬ 
vey (2MASS; Skrutskie et al. 1997). 2MASS characterize the 
large scale distribution of galaxies in the near-infrared using 
the K s (2.15/rm) passband. Our data set was selected from 
the public full-sky extended source catalogue (XSC; Jarret 
et al. 2000) which contains over 1.1 million extended objects 
brighter than K s = 14mag. The raw magnitudes where cor¬ 
rected for Galactic extinction using the IR reddening map of 
Schlegel, Finkbeiner & Davis (1998). There is a strong cor¬ 
relation between dust extinction and stellar density, which 
increases exponentially towards the galactic plane. Stellar 
density is a contaminant factor of the XSC since the relia¬ 
bility of separating stars from extended sources is very sen¬ 
sitive to this quantity (Jarret et al. 2000). In order to avoid 
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Figure 2. Redshift distribution of the different samples of clus¬ 
ters and galaxies analysed in this work. It can be seen that the 
distribution of the REFLEX clusters peaks at approximately the 
same redshift as 2MASS galaxies. 


contamination from stars we have constructed a mask for 
the 2MASS survey using a HEALPIX (Gorski et al. 1999) 
map with Nside = 256 and excluding those pixels where 
the K s band extinction A(K S ) = 0.367 x E(B — V) > 0.05 
and |6| > 20° which reduces galactic contaminant sources 
to 2% (Mailer et al., 2003). We also followed Mailer et al. 
(2003) and imposed a cut at K s = 13.57 mag in the cor¬ 
rected magnitudes. Figure 1 shows an Aitoff projection of 
the final sample obtained with these restrictions, which con¬ 
tains 397447 galaxies. 

We also used information from the APM Galaxy Sur¬ 
vey (Maddox et al. 1990), which is based on 185 UK Illa- 
J Schmidt photographic plates each corresponding to a 
5.8 x 5.8 deg 2 on the sky limited to bj ~ 20.5 and with 
a mean depth of ~ 400 /i _1 Mpc for b < —40° and 6 < —20°. 
We selected our galaxy sample from the APM survey im¬ 
posing the cut bj < 20.5. The sub-sample with bj < 18.5 is 
also shown in Figure 1. 

Galaxy K s luminosity is less sensitive to dust and stel¬ 
lar populations than the B band, providing a more uniform 
survey of the galaxy population. The K s band also gives a 
better measure of the stellar mass content than the B band. 
Then, the variations of the correlation function of samples 
selected in these bands can give us valuable information on 
the dependence of galaxy bias on the various properties of 
galaxies and then shed some light in their formation and 
evolution. 


2.2 Cluster and group catalogues 

In our analysis we used information of the ROSAT-ESO X- 
ray Flux Limited (REFLEX) Cluster Survey (Bohringer et 
al., 2004). The geometry of the survey is described by the 
southern hemisphere with S < +2.5° excluding the zone 
where |6| < 20° and the Small and the Large Magellanic 
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Clouds, covering an area of 4.24sr. The sample has a nom¬ 
inal flux limit of 3 x 10 -12 ergs -1 cm -2 within the ROSAT 
energy band (0.1-2.4)keV and comprises 447 clusters. Figure 
1 shows an Aitoff projection of the 422 clusters inside the 
2MASS mask described in section 2.1. 

We also used a group catalogue (hereafter the 2dF 
Galaxy Group Catalogue, 2dFGGC) constructed from the 
final version of the 2dF Galaxy Redshift Survey (2dFGRS, 
Colless et al. 2001), using the same technique used by 
Merchan & Zandivarez (2002) to construct the group cat¬ 
alogue of the 2dFGRS 100K release (Folkes et al. 1999). 
This sample comprises 5568 groups for which virial mass, 
velocity dispersion and virial radius have been determined. 

These samples where determined using very different 
types of information. The strong relation between Lx and 
M (see e.g. Reiprich & Bohringer, 2002) indicates that X- 
ray selected clusters samples are close to be basically se¬ 
lected by mass. The situation is very different for optically 
selected groups and clusters where the selection criteria is 
based on richness, which is not a good mass indicator. The 
study of the differences in the way in which galaxies cluster 
around them can be very important to improve our knowl¬ 
edge of this type of objects and to understand their overall 
statistical properties. 

Figure 2 shows the redshift distribution of the samples 
analysed in this work. The distribution for 2MASS galaxies 
was obtained by Mailer et al. (2004) by matching it with 
SDSS Early Data Release (Stoughton et al. 2002) and the 
2dFGRS 100K release (Colles et al. 2001). For APM galax¬ 
ies we used the fit from Baugh & Efstathiou (1993). It can 
be seen that the distribution for the REFLEX cluster sur¬ 
vey peaks at approximately the same redshift as 2MASS 
galaxies. As we measure the projected cross-correlation func¬ 
tion (see section 3), the largest contribution to this function 
will come for cluster-galaxy pairs laying at similar redshifts, 
which enhances our signal, on the other hand the differences 
with the distribution of APM galaxies will cause a lower am¬ 
plitude in the projected correlations. 


3 CLUSTER-GALAXY CORRELATIONS 

The spatial cluster-galaxy cross-correlation function £ cff (r) 
is defined so that the probability dP of finding a galaxy in 
the volume element dV at a distance r from the center of a 
cluster is 

dP = n[l+^ cg (r)]dV, (1) 


is the analogous quantity defined for a random distribution 
of Nr points covering the same angular mask than the orig¬ 
inal galaxy survey. 

When the angle 9 < 1 rad, that is, the distance y to 
the cluster centre is much larger than the projected cluster- 
galaxy separation, the projected correlation function is re¬ 
lated to the spatial correlation function by the simple in¬ 
tegral equation (Saunders, Rowan-Robinson & Lawrence, 
1992) 


£(0 = 


J_ f°° du cg (a) 1 
Bty J r da (a 2 — r 2 ) 1 / 2 


(3) 


The constant B regulates the amplitude of the correlation 
function taking into account the differences in the selection 
function of clusters and galaxies and can be calculated by 
(Lylje & Efstathiou 1988) 


B = 


EiV'(yi) 


E,N? fo°° S( x ) x2dx ’ 


(4) 


where ip(y) is the selection function of the galaxy survey, 
yt is the distance to cluster i and the sums extends to all 
clusters in the sample. 

Equation 3 can be easily solved analytically if we per¬ 
form a linear interpolation of uo cg between its values at the 
measured a’s. With this approach the solution to equation 
3 is (Saunders, Rowan-Robinson & Lawrence, 1992) 


£cg(ai) 


1 ©> +1 ~ ©> ln f °J +1 + V a J+i a i \ 

B7V i<j aj+1 ~ aj \ a i + ~ a i J 


We used this expression to obtain £ cg from ui cg . In order to 
calculate the factor B we evaluated the selection function 
S(v) by 

f- oo 

S(v) = / </>(L) dL, (6) 

J L min(y) 

where L m i n (y ) is the minimum luminosity that a galaxy at 
distance y must have in order to be included in the catalogue 
and <j>(x) is the luminosity function, for which we used the 
parameters of Bell et al. (2003) for 2MASS and Loveday et 
al. (1992) for APM. 


4 RESULTS 

4.1 The projected and 3D real-space 
cross-correlation function 


where h is the mean space density of galaxies. We only have 
information on the distance (redshifts) of the clusters and 
not on the individual galaxies in the 2MASS and APM cata¬ 
logues. Then, in order to obtain £ cg (r) we first determine the 
projected cross-correlation function ut cg (a), defined by Lilje 
& Efstathiou (1988), where a = cz6/Hq is the projected 
separation between a cluster at redshift z and a galaxy at 
angular distance 9 from its centre. To determine this quan¬ 
tity from our data we used the estimator 


u{a) 


(CG(a)) N r 
(CR( a)) No 


( 2 ) 


where Ng is the number of galaxies, ( CG(a )) is number of 
cluster-galaxy pairs at a projected distance a, and (CR(a)) 


The upper panels of Figure 3 show the projected cross¬ 
correlation function ui cg (a) obtained for REFLEX clusters 
and 2dF groups against APM (right) and 2MASS (left) 
galaxies for three different magnitude limits. In both cases 
the error bars where obtained by the bootstrap re-sampling 
technique. Specially for 2dF groups, the shape of this func¬ 
tion is different from a simple power law as it shows an in¬ 
crement in small scales that can be associated to the inner 
profile of the groups and clusters. 

The lower panels of Figure 3 show the corresponding 3D 
real-space cross-correlation function recovered from oj cg {a). 
Our findings are in complete agreement with the results of 
Croft et al. (1999, see their figure 2). We found that the 
shape of the correlation function can not be described as a 
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Figure 3. The upper panels show the projected cross-correlation function u> cg (cr) obtained for REFLEX clusters and 2dF groups against 
APM (right) and 2MASS (left) galaxies for three different magnitude limits. In both cases the error bars where obtained by the bootstrap 
re-sampling technique. The lower panels show the correspondent 3D real-space cross-correlation function recovered from uj cg (a). 


Table 1. Results for power-law fits of the form £ C g(r) = (r/ro ) 7 to the inner and outer regions of £ cg (r) for the different samples analysed 
using the REFLEX clusters as centres. 


REFLEX sample 

Galaxy sample 

pinner 

r 0 

pinner 

~outer 

r 0 

^outer 

all 

2MASS K s < 12.5 

6 1 +11 
° - 0.6 

—2 60~*~ 0 ' 26 
z.oU-o 17 

10 3+°- 90 
1U.O-1 17 

— 1 55 + 0,17 

1.OO-0 19 

all 

2MASS K s < 13.0 

°- z -0.7 

9 59 + 0-22 

9 85”*”°‘ 67 
y.oo_o 73 

— 1 69”*”°' 09 

i.oy _ 0 i4 

all 

2MASS K s < 13.57 

5 Q7 + 1 - 30 
-0.46 

—2 62” 1 ” 019 

9 is -1 " 0 - 60 
»-J-o_0.59 

i 77 +O.I 3 

A ** ' — 0.12 

all 

APM bj < 18.5 

5 47+0.64 
-0.38 

—2 69 +0 ' 24 
z,.u»_o io 


— 1 67+ 016 

1 - D *-0.15 

all 

APM bj < 19.5 

5 47 "*" 1,0 
-0.4 

— 2 61+0.25 
z - ol -0.17 

10.4±^f 

1 71+0-13 

± * /1 -0.14 

all 

APM bj < 20.5 

5 54 + 0-88 
0,0 -0.38 

—2 fi *)”*” 0 ' 21 

Z ’ DD -0.14 

9 5 +1-2 

— 1 79+°- 16 
<y -0.19 

log(Z,x/10 44 erg s _1 ) > —1.5 

2MASS Ks < 13.57 

5 96 +0 ' 79 
-0.46 

—2 62”*~ 0 ' 19 
Z - OZ -0.15 


1 77 +O.I 3 
-*-• 1 '0.13 

log(Lx/10 44 erg s _1 ) > —1.0 

2MASS K s < 13.57 

6 24 -1 " 0,95 
°- Z ^-0.61 

9 57+0-18 
-0.18 

io.28l° 0 ;g 

-1 70+°- 14 
' u —0.13 

log(Lx/10 44 erg s _1 ) > —0.85 

2MASS K s < 13.57 

6 24”*" 0-83 
O.Z4-0 57 

—2 [ 58”l” 0 ' 21 
Z.OO- 0.22 

10.45t°;®| 

1 75 + 0-10 
1 -‘°- 0.12 

log(Lx/10 44 erg s" 1 ) > —0.5 

2MASS K s < 13.57 

6 4i+°- 45 

u.^fci-o 27 

o 75 +O.H 
'°-0.06 

10 - 8 -0,7 

— 1 75 + 0,11 

A * *°-0.12 

log(Lx/ 10 44 erg s _1 ) > 0 

2MASS K s < 13.57 

7 1 q+ 0-89 
'■ lc> -0.44 

9 79 + 0.27 
' -0.07 


—1 69” 1 ” 015 

i.uy_ 0 is 


single power law for all scales, as it becomes steeper at small 
scales. Tables 1 and 2 show the results of power-law fits of 
the form £( r) cg = (r/ro) 7 to the inner and outer regions for 
the different samples analysed. It can be clearly seen that 
the correlation function is nearly independent of the magni¬ 


tude limit imposed to the galaxy sample showing only small 
variations in the best fit power-law parameters. The shape of 
the cross-correlations obtained for 2MASS and APM galax¬ 
ies are different. The feature at r ~ 2 /i -1 Mpc that show the 
transition to the inner profile of the clusters is more evident 
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Table 2. Results for power-law fits of the form £ cg (r) = (r/ro ) 7 to the inner and outer regions of £ cg (r) for the different samples analysed 
using 2dF groups as centres. 


2dF Groups sample 

Galaxy sample 

dinner 

r o 

pinner 

~outer 

r o 

ryOuter 

all 

2MASS K s < 12.5 

1 9 + 1 - 6 

iy -0.5 

0 90 +O .32 
°- ZO -0.26 

6-28j+j0 

— 1 60 + ° 07 
i.oU-o 07 

all 

2MASS K s < 13.0 

9 9 + 1.5 
Z ‘ Z -0.6 

9 70 + 0.34 
2.7o_o 35 

6 . 21 +°+ 

1 tr 1 + 0.06 
±,ol -0.06 

all 

2MASS K s < 13.57 

2 1+ 1 - 2 
Z-± -0.3 

—2 90”*” 0 ' 25 

0 40 

°' lo - 0.20 

1 tro + 0.06 
1 . 00-0 06 

all 

APM bj < 18.5 

I 0 + 1.1 
1 -®- 0.2 

1 7C+0-50 
± * '°-0.61 

6 fi ^ 4 " 9,33 

0 . 00-0 37 

1 07 + 0.08 
- L -°‘ -0.08 

all 

APM bj < 19.5 

2 6”*” 1 ' 3 
z>d -0.7 

9 r» 9 + 0.59 
-0.69 

c 71 +0.35 
°- / 1 -0.35 

— 1 42+ 0,09 
-*■•^-0.09 

all 

APM bj < 20.5 

Z-O - 0.6 

—2 20 +0 ' 43 

0.61 

4 82”*”°‘ 34 

0.36 

— 1 46”*” 0,12 
±, ^ D - 0.11 

log(A/„ ir /A/ 0 ) > 13 

2MASS K 3 < 13.57 

1 95+ 053 
0.16 

0 90 +O .82 
°' ZO -0.55 

6 73 "*” 0,27 

D -' -0.25 

1 [- 9 + 0.06 
1 ' OZ 0.07 

log(A/„ iT ./M Q ) > 13.5 

2MASS K s < 13.57 

9 9Q+0.65 

—3 47 + 0 - 94 
-0.31 

0 91 +0.39 
—0.35 

-1 59 +0 - 08 

-L.D»_o 07 

log(M„ ir /A/ 0 ) > 14 

2MASS K s < 13.57 

5-7 t\i 

—2 47 + 0 - 29 
-0.28 

10 . 6 ± 4 ° 

— 1 63” 1 ” 015 

1 .OO -0 17 

log(A/„ ir ./M 0 ) > 14.5 

2MASS K s < 13.57 

6-3+j 

—2 74 + 0 - 47 
'^-0.42 

10 *± 1:1 

— 1 84”*” 0,35 
l-o^_0.35 


for 2MASS galaxies. This effect can be associated with dif¬ 
ferences in the spatial distribution of these objects around 
clusters. 2MASS galaxies where selected in the near-infrared 
and so would tend to be luminous and of early type mor¬ 
phologies, and so, are expected to have a stronger clustering 
(Norberg et al. 2002). 


4.2 Dependence on Lx and M v i r 

In order to analyse the dependence of the distribution of 
galaxies around clusters and groups on properties such as 
their mass or X-ray luminosity we have computed £ cg (r) 
for 2MASS galaxies and the REFLEX cluster survey for 
four different limits on Lx- Analogously for 2dFGGC we 
have adopted four different limits in M v i r - The results ob¬ 
tained for these samples are shown in figures 4 and 5. We 
find a clear dependence of the amplitude of £ cg (r) on Lx 
which indicates that galaxies are more clustered around 
more luminous (massive) objects. The value of the ampli¬ 
tude of the correlaction function for the outer region changes 
from ro = 9.18 ± 0.59 for log(Lx/10 44 erg s -1 ) > —1.5 to 
ro = 13.1 ± 1.1 for log(Lx/10 44 erg s _1 ) > 0 and the inner 
region shows a similar behaviour. The shape of £ cg (r) does 
not show any important differences, but the small changes 
of ry lnner show a weak indication that a higher limiting Lx 
produce steeper inner regions. 

The changes in £ cg (r) for 2dF groups and 2MASS vary¬ 
ing the limiting virial mass are even stronger. The overall 
shape of £ , cg (r ) changes, specially for small scales where the 
amplitude increase from the sample with log {M v i r /M q) > 
13 to the one with log(A/„j r /M 0 ) > 14.5 is a factor 20 
whereas for scales r ~ 10 h^ 1 Mpc the increase is a fac¬ 
tor 5. The value of 7 in the outer region also increase with 
M v i r and for log (M„i r /Mg) > 14.5 its value is simmilar to 
the one in the inner region. 

Figure 5 also shows for comparison £ cg (r) for REFLEX 
clusters with log(Lx/10 44 erg s -1 ) > —0.85. Its similar¬ 
ity with the cross-correlation function obtained using 2dF 
groups with log (A/„j r /A/ 0 ) > 14 is remarkable. This fact 
can be used as an indirect test of the Lx — M relation. 
Reiprich & Bohringer (2002) found a power law relation be- 



j_1_1_ 

1 10 

r / (Mpc h _1 ) 


Figure 4. The real space cluster-galaxy cross-correlation func¬ 
tions of REFLEX clusters and all 2MASS galaxies for different 
X-ray luminosity limits. 


tween Lx in the ROSAT energy band (0.1-2.4 keV) and the 
mass A /200 (where A /200 = M tota i(r < r 2 00 )) of the form 


log 


_ Lx_ _ 

/ig O 2 10 4O erg s _1 


= A + a log 


f A /200 \ 
\ Ko ) 


(7) 


The transformation from Mvir to A /200 can be done assum¬ 
ing a given density profile for the clusters (see White 2000). 
For a ACDM the difference is small ( M v i r ~ 1.2A/2oo)- Then, 
assuming h = 0.7, A = —20.055 and a = 1.652 (from table 7 
in Reiprich and Bohringer, 2002), log(A/„j r //i , -1 Mq) = 14.0 
corresponds to log(Lx/10 44 erg s _1 ) = —0.85 in excellent 
agreement with the limit imposed to the sample used in fig¬ 
ure 5. 
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Figure 6. The results for b C g(r) for all the samples analysed in this work obtained by the comparison with the results of a ACDM N- 
body simulation. Different lines show the results obtained for samples selected with different limiting mass given by mi = 7.5 x 10 13 Mq 
( dot-dashed), m 2 = 1.4 x 10 14 M© (long dashed), m 3 =4.1 x 10 14 M© (short dashed) and 7714 = 7.5 x 10 14 M© (solid). 
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5 COMPARISON WITH THE RESULTS OF 
ACDM N-BODY SIMULATIONS 

In order to make a comparison between our results and 
those corresponding to the concordance cosmological model 
we used the ACDM Very Large Simulations (VLS). This 
simulations where carried out by the Virgo Supercomputing 
Consortium using computers based at Computing Centre of 
the Max-Planck Society in Garching and at the Edinburgh 
Parallel Computing Centre. The data are publicly avail¬ 
able at www.mpa-garching.mpg.de/NumCos. It uses 512 3 
particles in a box of 480 /U 4 Mpc, with initial conditions 
consistent with a ACDM power spectrum computed using 
CMBFAST (Seljak & Zaldarriaga, 1996), normalized so that 
as = 0.9. We have identified groups using a Friends-Of- 
Friends (FOF) algorithm with a linking parameter b = 0.2 
and we used these halos as centres to measure the halo-mass 
cross-correlation function £,hm(r) for different values of limit¬ 
ing mass. We used these results to infer the scale dependence 
of the cluster-galaxy bias factor, that is the bias between the 
distribution of galaxies and mass around clusters, defined by 

beg ( r ) = Cg ( r ) / £hm (r). (8) 

For small scales this is simply the difference between the 
density profile of mass and galaxies. 


In order to obtain b cg (r) we used four values of the lim¬ 
iting mass to measure £hm(r): mi = 7.5 x 10 13 Mq, m 2 = 
1.4 x 10 14 Mq, m 3 = 4.1 x 10 14 M 0 and m 4 = 7.5 x 10 14 Mq. 
As the REFLEX clusters and the 2dF groups span different 
ranges of mass, we used the 2dFGGC as centres to calculate 
Cg( r ) for m i and m 2 , and the REFLEX catalogue for m 3 
and m 4 . When using the 2dFGGC we simply imposed the 
same limits in m v i r . For the REFLEX clusters, we used the 
X-ray luminosity threshold for which the abundance, calcu¬ 
lated by the integration of the X-ray luminosity function of 
Bohringer et al. (2002), gives the same values as those ob¬ 
tained in the simulation. The results obtained are shown in 
Figure 6 . 

Our results for the two lower values of the threshold 
mass show that, at large scales, b cg (r) is consistent with a 
constant of order b cg ~ 1 for APM galaxies and slightly 
higher for 2MASS, which shows that 2MASS galaxies are 
more biased tracers of the mass distribution than optically 
selected galaxies. On smaller scales (r < 2 /i _ 1 Mpc) the 
distribution of galaxies around clusters shows a significant 
and nearly constant anti-bias of b cg ~ 0.7. For m3 and m.4, 
b cg ( r ) shows the same behavior but with a larger transition 
scale between the two regimes (r ~ 2 /U 4 Mpc). In all cases 
this scale would correspond to the onset of non-linearities, 
where £ cg (r) becomes steeper and virialization takes place. 
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Figure 5. The real space cluster-galaxy cross-correlation func¬ 
tions of 2dFGGC and all 2MASS galaxies for different lim¬ 
its on M v i r . For comparison £ cg (r) for REFLEX clusters with 
log(Lx/10 44 erg s _1 ) > —0.85 and 2MASS galaxies is also shown. 

A constant bias is a feature characteristic of linear theory 
and then it is not surprising that it does not hold as well in 
the non-linear regime. Another reason that may add to the 
anti-bias in the cluster and infall regions is the change in the 
galaxy types. The dense environments contain galaxies with 
older stellar populations with less light per baryonic mass. 
Then, cluster regions are down-weighted by this effect. It is 
very important to note that these results are independent of 
the limiting values of mass, X-ray luminosity and magnitude 
used to construct the clusters and galaxy samples, which 
shows that this is not a property of a given galaxy type, but 
a generic feature of the processes that control the efficiency 
of galaxy formation and evolution. 


6 DISCUSSION AND CONCLUSIONS 

In this work we have performed the first detailed calculation 
of the 3D real-space cluster-galaxy cross-correlation for an 
X-ray selected cluster sample and galaxies selected in the op¬ 
tical and near-infrared wave-bands. We have also presented 
the first calculation of this statistical measure for a galaxy 
group sample from the 2dFGRS. These two samples span a 
wide range of masses and so they allow us to analyse the 
dependence of £ cg (r) on M v i r . Our results for the different 
sub-samples of galaxies from 2MASS or APM with different 
limiting magnitudes show that £ cg (r) is almost independent 
of galaxy properties and that its shape is determined almost 
exclusively by the criteria used to define the cluster sample. 

In agreement with previous works (Croft et al. 1999) we 
found that the shape of £ cg (r) can not be described by a sin¬ 
gle power law at all scales. Instead, it shows two regimes with 
a clear transition, one for the larger scales, and a steeper one 
at smaller scales showing the inner profiles of the distribu¬ 
tion of galaxies around cluster centres. We have used our re¬ 


sults to check the Lx — M relation of Reiprich & Bohringer 
(2002). We have found that the correlation functions ob¬ 
tained using 2dF groups with a given lower mass limit and 
REFLEX clusters with the corresponding limit on Lx are 
in complete agreement showing the validity of this relation. 

The comparison of our results with those obtained for 
the halo-mass cross-correlation function in a ACDM N-body 
simulation shows that the observational results are consis¬ 
tent with a constant bias on large scales of order unity for 
APM galaxies and slightly higher for 2MASS galaxies. On 
smaller scales our results suggests that there is a substan¬ 
tial anti-bias ( b cg (r ) ~ 0.7). In all cases the transition scale 
between the two regimes corresponds to the onset of non- 
linearities, where the correlation functions becomes steeper. 
This is a strong result of our analysis which is independent of 
the properties of the cluster or galaxy samples used, showing 
that it is simply a generic feature of the different processes 
that govern galaxy formation and their subsequent evolu¬ 
tion. 
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